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Abstract 

We  apply  the  observability  rank  condition  to  the  two-dimensional 
Doi  orientation  tensor  model  to  investigate  the  observability  of  liquid 
crystalline  polymers  under  imposed  shear  or  extensional  flows.  Here 
the  observability  means  the  ability  of  determining  the  orientation  ten¬ 
sor  from  the  time  history  of  the  observations  of  the  birefringence.  The 
Doi  model  under  shear  or  extensional  flow  is  observable  almost  every¬ 
where.  In  other  words,  all  components  of  the  orientation  tensor  can  be 
determined  from  the  observation  of  the  birefringence  provided  that  the 
points  formed  by  the  two  components  of  the  orientation  tensor  do  not 
lie  on  certain  curves  in  the  plane.  We  also  run  simulations  using  the  un¬ 
scented  Kalman  filter  (UKF)  to  reconstruct  the  orientation  tensor  from 
observations  without  and  with  noises.  The  UKF  gives  good  estimates 
for  the  orientation  tensor  both  in  the  absence  and  in  the  presence  of 
observation  noises. 


1  Introduction 

Observability  of  a  dynamic  system  is  a  very  useful  subject  in  control  theory 
[3].  It  means  the  ability  to  determine  uniquely  the  state  of  the  system  from 
observable  quantities  [5,  8].  In  practice,  noise  may  occur  in  both  the  obser¬ 
vations  and  the  system  dynamics.  If  the  system  is  not  observable,  then  it  is 
hopeless  to  get  an  accurate  estimate  of  the  state  or  to  guide  the  system  to 
a  desirable  final  state.  So  from  a  practical  point  of  view,  it  is  important  to 
investigate  the  observability  of  a  dynamical  system. 

In  this  work  we  extend  our  previous  study  on  the  observability  of  viscoelas¬ 
tic  fluids  [8]  to  the  observability  of  the  Doi  orientation  tensor  model  in  liquid 
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crystalline  polymers.  We  impose  the  velocity  field  and  assume  monodomain 
structure.  We  further  assume  the  birefriengence  is  measured  in  experiments, 
which  is  consistent  with  light  scattering  measurements.  Our  goal  is  to  infer  the 
full  orientation  tensor  from  measurements,  if  possible.  To  do  so,  we  apply  the 
observability  rank  condition  to  investigate  the  short  time  local  observability  of 
the  Doi  orientation  tensor  model  under  imposed  extensional  or  shear  flow. 

The  outline  of  this  paper  is  as  follows.  We  present  a  brief  overview  of 
the  observability  of  dynamical  syterns  and  the  observability  rank  condition  in 
Section  2.  In  section  3  we  focus  on  the  analysis  of  the  observability  of  the 
Doi  orientation  tensor  model  under  shear  or  extensional  flow.  We  carry  out 
some  numerical  experiments  where  the  performance  of  the  unscented  Kalman 
filter  is  examined  both  in  the  case  of  zero  observation  noise  and  in  the  case  of 
Gaussian  observation  noise  in  Section  4.  We  summarize  our  work  in  Section 
5. 

2  Observability  rank  condition  of  dynamical 
systems 

For  reader’s  convenience,  we  quickly  recall  the  observability  of  dynamical 
systems  and  observability  rank  condition  [8]. 

Consider  a  dynamical  system  without  control  input: 


x  =  f(x),  x  G  Rn 

(1) 

y  =  Kx ),  y  ^RP 

(2) 

a;(0)  =  a:0, 

(3) 

where  x  is  the  state  variable,  y  is  the  quantity  that  can  be  observed  experi¬ 
mentally.  Functions  /  and  h  are  assumed  to  be  sufficiently  smooth.  The  state 
x  is  not  observed  directly  but  the  output  y  is.  We  would  like  to  know  whether 
it  is  possible  to  determine  the  initial  state  x°  from  the  output  history  y{ 0  :  oo) 
where  the  symbol  y( 0  :  T)  denotes  the  trajectory  t  — >  y(t)  with  0  <  t  <  T. 
If  the  map  a;0  — »  y( 0  :  oo)  is  one  to  one,  then  one  can  reconstruct  x°  from 
y( 0  :  oo)  mathemematically  and  the  system  (l)-(2)  is  called  observable.  Since 
x{t)  can  be  treated  as  the  initial  state  for  evolutions  after  time  t,  the  observ¬ 
ability  indicates  that  x(t )  can  be  reconstructed  from  y(0  :  oo)  as  well. 

Other  definitions  of  observability  are  also  possible.  The  most  useful  and 
easiest  to  measure  one  is  short  time  local  observability.  A  system  is  short  time 
locally  observable  if  for  every  T  >  0,  the  map  x°  — >  y( 0  :  T )  is  locally  one 
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to  one.  Mathematically  speaking,  locally  one  to  one  means  that  for  each  x°, 
there  is  a  neighborhood  U(x°)  such  that  if  a:1  G  U(x°),  then  the  output  from 
xl  and  the  output  from  x°  are  different.  Now  we  recall  a  sufficient  condition 
for  the  short-time  local  observability  [5]. 

We  begin  with  some  notations.  Note  that  from  (2)  and  (3)  we  have  y(0)  = 
h(x°)  and  thus 

f  (0)  =  <4> 

using  the  chain  rule  and  (1).  The  right  hand  side  of  (4)  is  called  the  Lie 
derivative  of  the  function  h  by  the  vector  held  /: 

dh 

Lf{h){x)  =  —{x)f(x).  (5) 

Since  the  Lie  derivative  is  another  function  from  Rn  to  Rp,  one  can  repeat  the 
process  of  taking  the  Lie  derivative: 

f)h 

L) )  (h)(x)  =  Lf(Lf{h))(x)  =  LfQ(x)f(x)) 

Lk(h)(x)  =  Lf(Lkf\h))(x). 

A  set  of  functions  gi(x),  ■  ■  •,  gk(x)  is  said  to  separate  points  if  given  any  pair 
of  two  points  x°  and  x1,  there  is  at  least  one  gi(x)  such  that  gi(x°)  ^  gi{xl).  If 
g i(x),  ■  ■  •,  gk(x)  separate  points,  then  the  map  x  — >  (gi(x),  ■  ■  ■ ,  guix))  is  one  to 
one.  Mathematically,  if  the  functions  Lj.  (h)(x),  j  —  0,  •  •  • ,  k,  locally  spearate 
points  in  Rn  for  some  k ,  then  the  system  is  short  time  locally  observable.  A 
sufficient  condition  for  this  to  happen  is  that  the  so-called  one  forms 

dhi(x),  •  •  • ,  dLk  (hi)(x),  i  =  !,-•*, p 


span  n  dimensions  at  every  x  where 


Now  we  introduce  the  definition  of  observability  rank  condition  (ORC). 

Definition  2.1  The  system  (l)-(2)  satsifies  the  observability  rank  condition  at 
x°  if  there  exists  a  k  such  that  { d,L (hf)  :  j  —  0,  •  •  k\  i  —  1,  •  •  •  ,p}  has  rank 
n.  The  system  (l)-(2)  satisfies  the  observability  rank  condition  if  it  satisfies  it 
at  every  x  G  Rn  (note  k  may  vary  with  x). 
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For  linear  systems,  observability  rank  condition  (ORC)  leads  to  global 
short-time  observability  as  well.  For  nonlinear  systems,  ORC  is  a  sufficient 
condition  of  short-time  local  observability.  In  addition,  ORC  is  almost  a 
necessary  condition  for  short-time  local  observability.  It  has  been  shown 
HeTm&rmKrennerthatiftheORCisviolatedonanopensubsetofW1 ,  then  the  sys¬ 
tem  (l)-(2)  is  not  short  time,  locally  observable. 

In  the  next  section  we  would  like  to  apply  the  observability  rank  condition 
to  investigate  the  observability  of  the  Doi  orientation  tensor  model  in  the  study 
of  liquid  crystalline  polymers. 


3  The  Doi  orientation  tensor  model 

The  Doi  orientation  tensor  model  for  rodlike  liquid  crystalline  polymers  in 
a  solvent  is  well-known  for  its  capability  to  describe  both  the  isotropic  and 
nematic  phases  and  phase  transition  between  them  [6].  A  fundamental  part  of 
this  model  is  the  single  molecule  orientation  distribution  function.  Interactions 
between  molecules  are  represented  by  a  mean-held  potential  (so-called  Maier- 
Saupe  potential).  The  rodlike  molecules  are  also  subject  to  Brownian  force 
due  to  the  fact  that  they  interact  with  other  rodlike  molecules  and  with  the 
how. 

Basically,  the  model  is  a  microscopic  Smoluchowski  equation  or  Fokker- 
Planck  type  equation  for  the  dynamics  of  the  orientational  distribution  func¬ 
tion  coupled  with  a  macroscopic  hydrodynamic  equation  [1],  The  Smolu¬ 
chowski  equation  describes  the  convection,  rotation  and  diffusion  of  the  rodlike 
molecules.  The  full  Doi  orientation  tensor  theory  is  developed  after  the  kinetic 
Smoluchowski  equation  is  projected  onto  a  second-moment  description  using 
various  closure  rules.  The  major  element  in  this  tensor  theory  is  the  second- 
moment  tensor  Q  which  describes  the  orientational  distribution  of  the  ensemble 
of  rodlike  macromolecules.  The  orientation  tensor  is  traceless  and  symmetric. 
The  physical  and  practical  signihcance  of  the  orientation  tensor  is  that  it  is 
the  basis  for  micro-scale  light  scattering  measurements  of  primary  axes  (“di¬ 
rectors”),  degrees  of  molecular  alignment  (“birefringence”),  and  normal  and 
shear  stress  measurements. 

Here  we  restrict  our  attention  to  two-dimensional  monodomain  structure. 
The  study  of  two-dimensional  liquid  crystal  polymers  has  been  physically  in¬ 
spired  by  monolayer  films.  Thin  films  of  liquid  crystal  polymers  have  been 
used  as  alignment  layers  for  liquid  crystal  displays  because  of  their  stability 
and  nonlinear  optical  properties.  A  lot  of  theoretical  and  experimental  studies 
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have  been  devoted  to  the  two-dimensional  Doi  model  (for  example,  see  [7]  and 
references  therein). 

The  two-dimensional  Doi  model  is  given  by  [7] 

—  Q  =  12Q  —  Qf2  +  a[DQ  +  QD]  +  oD  —  2aD  :  Q(Q  +  ~)  —  6-DrU(Q),  (7) 

LLL  z 

where  Q  is  the  2x2  orientation  tensor  such  that 

Q  =  ^12  )  Qu  =  —Q22,  Q 12  =  Q21, 

V21  V22 

n  is  the  vorticity  tensor  such  that  12  =  |[Vv  —  VvT];  a  is  a  dimensionless 
parameter  which  depends  on  the  molecular  aspect  ratio;  D  is  the  rate-of-strain 
tensor  D  =  |[Vv  +  VvT];  Dr  is  the  averaged  rotary  diffusivity  or  relaxation 
rate;  F(Q)  is  defined  by 

F( Q)  =  (1  -  j)Q  -  NQ2  +  JVQ  :  Q(Q  +  1) 

and  N  is  a  dimensionless  concentration  of  nematic  polymers. 


3.1  The  Doi  model  under  extensional  flow 

When  an  extensional  flow  with  rate  7  is  imposed  to  the  monodomain  of 
liquid  crystalline  polymers,  the  velocity  field  can  be  expressed  as 

v  =  (7(f)f-.-7(f)|)»  (8) 

and  the  velocity  gradient  is 

7  (*) 

2 

0 

After  applying  the  gradient  of  the  homogeneous  extensional  flow  (9)  to  the  Doi 
system  (7),  we  obtain 


Qu  —  —6 Dr  [(1  —  y)Q  11  +  2NQ11(Qf1  +  Qf2)]  +  a(|  —  2Qli)^{t), 
Q12  —  —6 DrQi2  [(1  —  y)  +  2N(Qi1  +  Q12)]  ~  2aQuQi2'y(t). 


Using  the  nematic  relaxation  time  scale  A-,  the  flow  field  and  orientation  dy¬ 
namics  of  (7)  can  be  non-dimensionalized.  The  key  dimensionless  parameters 
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are  then  the  Peclet  number  Pe  =  7 / Dr  (the  shear  rate  normalized  with  respect 
to  nematic  relaxation  rate)  and  the  dimensionless  concentration  parameter  N. 
Rescaling  time  as  t  =  tDr ,  the  nematodynamic  model  (7)  in  the  dimensionless 
form  becomes 


Q11  —  —QQ11  [l  —  y  +  2  N(Q‘f1  +  Q12)]  +  Q  ( |  —  2Qf1)Pe, 
Q 12  =  —6Q12  [l  —  y  +  27V(Qi1  +  Qf2)]  —  2aQuQi2Pe. 


More  succinctly,  the  system  (11)  can  be  written  as 


where 


X\ 

Q  ii 

.  x2  . 

Q 12 

—6  1  —  ^  +  2N{x\  +  xf)  x  +  aPe 


2xj 


-2x\x2 


(12) 


In  the  light  scattering  experiments  one  can  measure  the  birefringence,  which 
corresponds  to  the  difference  of  the  two  eigenvalues  of  Q.  Since  the  two  eigen¬ 
values  of  Q  are  Ai;2  =  ±4/ Q h  +  Q\2,  their  difference  is  2a/ Qh  +  Q\2.  For 
computational  convenience,  we  assume  the  observation  is  h  =  Q'\{  +  Q\2  = 

ry»  2  I  rp  2 
<X/  |  <X/  2  • 

The  Lie  bracket  is 


dh_ 

dx.i 


(x)fi{x) 


"  dh  dh ' 

'  /i  ' 

dxi  ’  (9x2 

.  /2  . 

=  [2x1,2x2 


h 

/2 


—  2  (xi  fi  +  x2  f2) 

=  —  12x\  +  61Vxf  —  241Vxf  —  4&Nx\x\  +  aPex  1 
—4  aPex\  —  I2x2  +  QNx\  —  2ANx\  —  4aPexix| 


After  some  algebra,  one  has 


f  - 

h 

^  =  det 

2xi 

2x2 

.  Lf^  . 

ci 

C2 

=  2aPex2(Ax\  +  4x2  —  1), 


(13) 


(14) 
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where 

Ci  =  —24^1  +  l2Nxi  —  96Nxl  —  9QNx\x\  +  a  Pe  —  12 a  Pex\  —  4 a  Pex |, 
c2  =  —9§Nx\x2  —  24x2  +  12Nx2  —  9QNx\  —  8aPexix2- 

(15) 

So  if  x\  +  x%  7^  1/4  and  x2  7^  0,  then  the  determinant  is  nonsingular  and 
consequently  the  observability  rank  condition  (ORC)  is  satished.  In  other 
words,  if  the  two  components  of  the  orientation  tensor  Q  satisfies  the  conditions 
Qi  i  +  Q\ 2  7^  1/4  and  Q\2  7^  0,  then  the  Doi  model  under  extensional  flow  is 
short  time  locally  observable. 

In  Figure  1  we  plot  the  two  curves  described  by  x2  =  0  and  x\  +  x\  =  1/4, 
respectively  in  the  X\  —  x2  plane..  The  Doi  model  under  extensional  flow  is 
short  time,  locally  observable  as  long  as  the  point  {Qii,Q\2)  does  not  lie  on 
these  two  dashed  curves. 
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Figure  1:  The  curves  where  the  determinant  in  (14)  vanishes.  The  Doi  model 
under  extensional  flow  is  observable  for  all  components  Q11  and  Q\2  of  the 
orientation  tensor  when  the  point  (QiuQ  12)  doesn’t  lie  on  these  curves. 
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3.2  The  Doi  model  under  shear  flow 

The  velocity  field  of  a  shear  flow  with  rate  7  is  described  by 

v  =  i(t)(y,0) 

The  rate-of-strain  tensor  is 

'  0  1 


D  =  i[Vv  +  (Vv)r]  = 


1  0 


(16) 


(17) 


We  non-dimensionalize  the  flow  field  and  orientation  dynamics  using  the  ne¬ 
matic  relaxation  time  scale  1  / Dr.  The  dimensionless  parameters  consist  of  the 
Peclet  number  Pe  =  7 (t)/Dr  which  is  the  shear  rate  normalized  with  respect 
to  nematic  relaxation  rate,  and  the  dimensionless  concentration  parameter  N. 
From  now  on,  we  work  in  dimensionless  time  t  =  tDr. 

The  nematodynamic  Doi  model  (7)  in  component  form  is 

Qu  =  —QQ11  [l  —  y  +  2-/V(Qi!  +  Q\ 2)]  +  (Q12  —  2aQuQi2)Pe, 

Q 12  =  —6Q12  [l  —  y  +  2,N(Qf1  +  Q12)]  +  (— Qu  +  f  —  2aQ\2)P e. 
Introducing 


(18) 


x  = 


X\ 

Q 11 

.  X2  . 

Q 12 

the  system  (18)  has  the  simple  form 

dx 


dt  =  /(f)’ 

f(x)  = 


h 

/2 


=  -6  [1  -  y  +  2N(x\  +  x|)]  x  +  Pe 


x2  —  2axix2 


-xi  +  f 


2 ax2  j 


(19) 

As  before,  assume  the  the  observation  is  related  to  the  birefringence  and 
has  the  form  h  =  Q2U  +  Q\2  =  x\  +  x\  . 

The  Lie  bracket  is 

(*)/»(*)  = 


'dh  dh' 

'  fi  ' 

dxi  ’  cLc2 

.  h  . 

=  [2xi,2x2 


'  fi  ' 

.  h  . 


2  (xi  fi  +  x2  f2) 


=  —  12x\  +  6iVxf  —  24  Nx\  —  4%Nx\x1  —  4aPex\x2 


— 12^2  +  &Nx2  —  24Nx2  +  x2  a  Pe  —  4 a  Pex\. 


(20) 
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It  follows  that 


det  | 

ox 


h 

Lfh 


=  det 


2x\  2x2 
ch  d  2 


=  —  2aPexi(4:Xi  +  Ax2  —  1),  (21) 


where 

di  =  —2Axi  +  12Nx±  —  96Nxl  —  96Nxyxl  —  8aPexix2, 


d2  =  —9bNx\x2  —  4  aPex\  —  24x2  +  12Ahc2  —  96Nx%  +  aPe  —  12  aPex^. 

(22) 

Therefore,  if  X\  ^  0  and  x\  +x\  ^  1/4,  then  the  matrix  is  nonsingular  and  the 
ORC  is  satisfied.  Said  differently,  if  Q n  ^  0  and  Qh  +  Q\2  ^  1/4,  then  the 
Doi  model  under  shear  flow  is  short  time  locally  observable. 

We  plot  the  two  curves  described  by  x\  =  0  and  x\+x\  =  1/4,  respectively 
in  the  X\  —  x2  plane  in  Figure  2.  The  Doi  model  under  shear  flow  is  short 
time,  locally  observable  as  long  as  the  point  (QmQ  12)  does  not  lie  on  these 
two  curves. 
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Figure  2:  The  curves  where  the  determinant  in  (21)  becomes  zero.  The  Doi 
model  under  shear  flow  is  observable  for  all  components  Q n  and  Q i2  of  the 
orientation  tensor  when  the  point  (QiuQ  12)  doesn’t  lie  on  these  curves. 
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4  Unscented  Kalman  filtering  (UKF)  of  the 
Doi  model 


The  observability  of  a  dynamic  system,  as  studied  above,  provides  the 
mathematical  possibility  of  recovering  the  state  of  the  system  from  observable 
quantity.  However,  the  concept  of  observability  itself  does  not  provide  the 
mechanism  to  practically  extract  the  state  from  observations.  In  contrast, 
Kalman  filter  is  a  widely  used  tool  to  estimate  the  state  using  observation 
data  for  systems  with  observability. 

Now  we  want  to  investigate  the  performance  of  the  unscented  Kalman  fil¬ 
tering  [4]  in  recovering  the  orientation  tensor  from  the  observed  birefringences. 
We  refer  the  readers  to  the  original  paper  by  Julier,  Uhlmann  and  Durrant- 
Whyte  [4]  for  a  complete  description  of  the  unscented  Kalman  filtering.  UKF 
is  a  popular  choice  when  the  state  transition  and  observation  models  are  highly 
nonlinear.  However,  a  rigorous  proof  on  its  convergence  is  still  open.  In  this 
section  we  will  examine  two  cases:  (1)  the  observation  is  noiseless  and  (2)  the 
observation  is  polluted  by  Gaussian  noise.  We  select  the  Doi  model  under 
shear  flow  as  the  testbed.  Similar  results  are  obtained  for  the  Doi  model  under 
extensional  flow  which  are  omitted  here. 

To  examine  UKF  for  a  non-steady  state  solution  of  the  Doi  model  under 
shear  flow,  we  add  an  external  force 


0.1  sin  (t) 
0.1  cos(t) 


(23) 


to  the  right-hand  side  of  equation  (11).  We  remark  that  the  observability  study 
in  the  previous  section  is  still  valid  here. 

We  start  the  system  of  Doi  model  under  extensional  flow  (11)  with  the 
external  force  (23)  at  some  initial  condition  T(0)  and  the  filter  at  a  different 
initial  condition  x(0).  We  solve  the  system  equations  without  noise  to  get  state 
T(0  :  ex))  and  observation  h( 0  :  oo)  trajectories.  In  the  case  of  observations 
with  no  noise,  we  pass  the  noise  free  observation  trajectory  to  the  filter  and 
the  filter  yields  a  state  estimate  trajectory  x(()  :  oo).  The  estimation  error 
x(t)  —  x{t)  is  the  difference  between  the  state  of  the  system  (which  is  not 
directly  measurable)  and  the  estimate  state  produced  by  the  filter  from  the 
observation.  The  filter  is  said  to  be  convergent  if  the  estimation  error  goes  to 
zero  as  t  — >  oo  for  any  T(0)  and  £(0). 

In  Figure  3  we  show  the  results  of  UKF  for  the  Doi  model  under  shear  flow. 
The  parameters  used  here  are  N  =  4.0,  Pe  =  0.01,  and  a  =  0.8.  Figure  3(a) 
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depicts  the  exact  and  estimated  solutions  of  the  Doi  system  under  shear  flow, 
respectively.  Figure  3(b)  gives  the  corresponding  filter  errors.  It  implies  the 
convergence  of  the  UKF  when  there  is  no  noise  but  only  initial  estimate  error. 


(a) 


(b) 


Figure  3:  (a)  Exact  solutions  (solid  lines)  and  estimated  solutions  from  UKF 
(dashed  or  dotted  line)  of  the  Doi  model  under  shear  flow  with  an  external 
force,  (b)  Corresponding  filter  errors  of  UKF  in  (a). 


In  practice  all  observations  may  be  polluted  by  noises  from  different  sources. 
Next  we  check  the  performance  of  UKF  in  the  presence  of  observation  noise. 
We  add  Gaussian  noise  to  the  measurements.  Specifically,  we  assume  that  the 
experimental  measurement  is  the  true  value  of  the  observation  plus  Gaussian 
noise: 


^experiment  (0  :  °°)  =  MO  :  °°)  +  Gaussian  noise. 

Then  we  pass  the  observation  with  noise  ^experiment  (0  •  °°)  to  the  filter. 
The  filter  outputs  an  estimated  state.  Figure  4  depicts  the  results  of  UKF 
for  the  Doi  model  under  shear  flow.  Figure  4(a)  gives  the  measurement  with 
noises  (the  solid  line)  vs  the  measurement  without  noise  (dashed  line)  whereas 
Figure  4(b)  isolates  the  noise  added  to  the  observation.  The  exact  solutions  and 
estimated  solutions  from  UKF  are  plotted  in  Figure  4(c)  and  the  corresponding 
filter  errors  are  given  in  Figure  4(d).  It  is  clear  that  UKF  yields  good  estimate 
even  when  noises  are  present  in  the  measurements. 
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5  Conclusions 

We  have  applied  the  observability  rank  condition  to  the  vector  holds  in  the 
Doi  orientation  tensor  model  to  investigate  the  short  time  local  observability 
of  liquid  crytalline  polymers  driven  by  extensional  or  shear  how  helds.  The 
measurement  is  assumed  to  be  birefringence  of  the  materials  which  is  consistent 
with  experiments.  We  have  found  that  the  Doi  model  is  observable  almost 
everywhere.  That  is,  one  can  determine  the  orientation  tensor  from  the  time 
history  of  the  observations  of  birefringence  when  the  points  formed  by  the  two 
components  of  the  orientation  tensor  do  not  lie  on  certain  curves  in  the  plane. 
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(C)  (d) 


Figure  4:  (a)  Observation  with  Gaussian  noise  (solid  line)  vs  observation  with¬ 
out  noise  (dashed  line),  (b)  The  Gaussian  noise  in  (a),  (c)  Exact  solutions 
(solid  lines)  and  estimated  solutions  from  UKF  (symbols)  of  the  Doi  model 
under  shear  flow  with  an  external  force  where  observations  contain  Gaussian 
noise,  (d)  The  corresponding  filter  errors  of  (c). 


